Compressible dynamics of magnetic field lines for 
incompressible MHD flows 

E.A. Kuznetsov^t T. Passot ^ and P.L. Sulem W 

( a ) - Landau Institute for Theoretical Physics, 2 Kosygin str., 119334 Moscow, Russia 
W - CNRS, Observatoire de la Cote d'Azur, PB 4229, 06304 Nice Cedex 4, France 



Abstract 

It is demonstrated that the deformation of magnetic field lines in incom- 
pressible magnetohydrodynamic flows results from a compressible mapping. 
Appearance of zeroes for the mapping Jacobian correspond to the breaking of 
magnetic field lines, associated with local blowup of the magnetic field. The 
possibility of such events is found to be unlikely in two dimensions but not in 
three dimensions. 

PACS: 52.30. Cv, 47.65.+a, 52.35.Ra 



1 Introduction 

An important property of ideal magnetohydrodynamics (MHD) is the frozenness of 
magnetic field in the plasma: fluid particles remain pasted on their magnetic lines 
that are driven by the transverse velocity component. This remark is the start- 
ing point of a mixed Lagrangian-Eulerian description of ideal MHD flows, named 
magnetic line representation (MLR) and first formulated in [1] . The idea originates 
from the vortex line representation (VLR) of hydrodynamic flows [2] that involves 
a two-dimensional Lagrangian marker labeling each vortex line, together with a pa- 
rameterization of these lines. In three dimensions (3D), this representation enables 
one to partially integrate the Euler equations with respect to a continuous infinity of 
integrals of motion called the Cauchy invariants. A main peculiarity of the transfor- 
mation associated with the vortex line dynamics is its compressible character that, 
as recently pointed out by one of the authors [3], is amenable of a simple interpreta- 
tion. The Euler equations can be rewritten as the equations of motion for a charged 
compressible fluid moving under the action of effective self-consistent electric and 
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magnetic fields satisfying Maxwell equations. The new velocity coincides with the 
velocity component transverse to vorticity, which, due to the frozenness property, 
identifies with the vortex line velocity. It is well known that the appearance of 
singularities in compressible flows is connected with the emergence of shocks, cor- 
responding to the formation of folds in the classical catastrophe theory [4]. In the 
gas-dynamic case, this process is completely characterized by the mapping defined 
by the transition from the usual Eulerian to the Lagrangian description. A zero of 
the Jacobian corresponds to the emergence of a singularity for the spatial deriva- 
tives of the velocity and density of the fluid. Due to the compressible character of 
VLR, the phenomenon of breaking becomes also possible for vortex lines in ideal 
incompressible fluids. Vortex line breaking was first studied for three-dimensional 
integrable hydrodynamics with Hamiltonian Ji — j \ Q\dr where ft is the vorticity 
[5]. This model and the Euler equation are both incompressible and have the same 
symplectic operator defining the Poisson structure. Breaking of vortex lines is as- 
sociated with the touching of two vortex lines and results in an infinite vorticity. 
Recent numerical simulations [6, 7] have suggested the possibility of such a scenario 
for the 3D Euler equations, but further investigations are required to reach a defi- 
nite conclusion. In ideal MHD, we can expect the same behavior for the magnetic 
field which is a frozen-in quantity. In two dimensions (2D) however, the fact that 
vorticity is perpendicular to the flow plane while the magnetic field lies in it, puts 
a limit to the analogy, making magnetic field line breaking a priori possible in two 
dimensions, while singularities are excluded in 2D Euler flows. It will nevertheless 
be argued in this paper that magnetic field blowup is unlikely in 2D MHD. 

In Section 2, we recall the Cauchy formula for MHD flows, which plays a central 
role in the derivation of the Weber type transformation discussed in Section 3. This 
transformation is obtained by extending ideas of paper [3] to ideal incompressible 
MHD flows. We in particular indicate how the MHD equations can be partially 
integrated. Section 4 addresses the two-dimensional case where two conservation 
laws are established. In Section 5, we discuss the possibility of magnetic line break- 
ing as a local blowup of the magnetic field, a process different from the gradient 
singularity associated with current sheets formation ([8] and references therein). A 
brief conclusion is provided by Section 6. 

2 Cauchy formula in MHD 

As well known, the magnetic field h in ideal incompressible MHD obeys 



that formally coincides with the equation governing the vorticity f2 in Euler hy- 
drodynamics. Since only the transverse velocity vj_ to the local magnetic field is 
relevant in this equation, we introduce new Lagrangian trajectories 



h t = curl(v x h), divv = 0, 



(1) 




(2) 
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defined by 

!=v l( r, t ) (3) 
r|t=o = a. (4) 

It is easily established that the Jacobian matrix (of element J^- = -^-) obeys 

^J = JU (5) 



where the matrix U has elements Uij = - 3 . One then obtains the equations for 



dv±j 

the Jacobian J = det J and for the inverse matrix J -1 with elements dctj/dxi, 
(where a = a(r, t) is the inverse of the mapping defined in (2)), in the form 



— J = div v_|_ J (6) 

(JjL 

and 

^.j-i = -uj-\ (7) 

Since divv^ is generically non zero, the mapping (2) is compressible and the Jaco- 
bian J can vanish. This observation is central in the discussion of the possibility of 
magnetic field blowup presented in Section 5. 

By means of eqs. (6) and (7), eq. (1) is transformed into 

a = 0, (8) 

where D t = d t + (vj_ • V) identifies with the material derivative d/dt used in (3). 
Integration of this equation leads to a "new" vector Lagrangian invariant 

/,(a) = J^ (9) 

that coincides with the initial magnetic field h (a) and is the analog of the Cauchy 
invariants for ideal hydrodynamics. The magnetic field h is then given by 

h(MH (h (a) • V a )r(a,Q (1Q) 

J 

3 Weber type transformation 

Equation (10) is the basis of the magnetic line representation [1]. Another important 
formula for MLR follows from the velocity equation 

d t w + (v • V)v = -Vp + curl h x h, (11) 
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where we normalized the magnetic field by the factor \/Airp (where p is the density) 
so that h has the dimension of a velocity 

We also decompose the velocity v = vj_ + v r into transverse and tangential 
components and substitute in (11). As a result, eq. (11) is rewritten as 

<9fV_|_ + (v_|_V)v_i_ = E + v_|_ x H + F mhd , (12) 

where we introduced new effective "electric" and "magnetic" fields 

(13) 

H = rot v T . (14) 
In eq. (12), the force F mhd = j x h, involves the renormalized current 

j = curl h — (v • h)/h 2 curl v. (15) 

As seen from (13) and (14), the new auxiliary "electric" and "magnetic" fields 

v 2 

can be expressed in terms of scalar and vector potentials •p = p + — and A = v r , 
so that the two Maxwell equations 

div H = 0, — = -curl E 

at 

are automatically satisfied. In this case, the vector potential A has the gauge 

div A = —div v_l, 

which is equivalent to the incompressibility condition div v = 0. 

The two other Maxwell equations define auxiliary charge density and current 
which follow from relations (13) and (14). 

New terms in the right hand side of eq. (12) also have a mechanical interpreta- 
tion. The Lorentz force vixH plays the role of a Coriolis force. The potential <p 
has a direct connection with the Bernoulli formula. The term d t v T results from the 
non-inertial character of the coordinate system. 

In eq. (12), we make the change of variable defined by mapping (2). As a result, 
the equations of motion are expressed in a quasi-Hamiltonian form, analogous to 
eq. (20) of [3] 1 

AP = -f + F^, Ar = |£, (16) 
or oP 

where the Hamiltonian h is given by the standard expression 

/ i = i(P-A) 2 + ^ = p+y, 



1 The first equation of the system (16) contains an addition term F mhd and therefore we qualify 
(16) of quasi-Hamiltonian. 
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in terms of the generalized momentum P = vj_ + A (that identifies with v), and 
thus coincides with the Bernoulli "invariant" for a non-magnetic fluid. 
Introducing a new vector 

Uk = Pi 7^—, 

depending on t and a, one easily obtains from (16) that this vector obeys 



Using (10) and the identity 

one has 
where 

Equation (17) thus rewrites 



dxidxj da a 
dap oa y ox k 



F rnM^ = h()(a) x g) 



S = (j • V r )a. 



D t u = V a (^-p+ V f- V -pj+ h (a) x S. (19) 

Integrating in time then leads to the Weber type transformation 

u = u (a) + V a Hh (a)xW, (20) 
where the potential $ satisfies a Bernoulli type equation, 

F 2 2 

and the vector W obeys 

D t W = S. (21) 

If initially $|t=o = and W| t=0 = 0, the integration "constant" u (a) coincides 
with the initial velocity v (a). This vector u (a) is thus a new Lagrangian invariant. 

To get a closed description we eliminate the pressure p by applying the curl 
operator (with respect to a-variables) on eq. (20) 

curl a u = curl a u (a) + curl a [ho (a) x W]. (22) 

This equation can also be rewritten as 

n(r,t) = (n ° (a ' t) - V ° )r(a ' t) . (23) 
u 
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Here fi (a,£) is given by 

n (a,t) = fio(a) + curl a [ho (a) x W], 

where f2o(a) is the initial vorticity. When ho (a) = 0, eq. (23) reduces to the Cauchy 
formula for vorticity in ideal hydrodynamics. 

The vector W is determined from eq. (21) that rewrites 

AW =^H,(M)-jcurl.(^) (24) 
where g is the MLR metric tensor defined by 

^ a/3 da a dap 

and b = Jh is given by (10). 

As a result, we have two equations of motion for the mapping (3) and for the 
vector W. Together with eqs. (10), (23) and the relation between velocity and 
vorticity, 

ft = curl r v, div r v = 0, (25) 

this constitutes a closed system of equations that provides a magnetic line repre- 
sentation for incompressible MHD (to be compared with [1]). These equations are 
solved with respect to two Lagrangian invariants h (a) and f2o( a )- It is possible to 
show [1] that conservation of these invariants in MHD is a consequence of relabeling 
symmetry, as it is the case for Euler equation (see, e.g. the reviews [9, 10]). 

The magnetic line representation involving the local change of variables r = 
r(a, t), breaks down at singular points where the Jacobian is zero or infinity and 
the normal velocity is not defined. 

Let us consider the null point r = r(t) defined by 

h(r(t),i) = 0. (26) 
Differentiating this equation with respect to time, we get 

^ + (r(t) • V)h = 0, 

with r(t) = v(r(t),t), which shows that the null points are advected by the flow. 
The velocity v at these points is defined by inverting the curl operator in (25). 

Null-points are topological singularities for the tangent vector field r(r). Their 
classification depends on the space dimension D. Topological constraints that can be 
considered as additional conditions for the MLR system, can be written as integrals 
of the vector field r(r) and its derivatives over the boundary of simply-connected 
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regions (in 3D) or along closed contours (in 2D) enclosing the null-points. In D = 2, 
one has 

£(Vip ■ dr) = 2nm, (27) 

where ip is the angle between the vector r and the x-axis and m is an integer often 
called topological charge. It is equal to the total number of turns of the vector r 
while passing around the closed contour encircling the null-point. For instance, for 
X-points or O-points, m = ±1. 

In D = 3, the topological charge is defined as the degree of the mapping S 2 — > S 2 , 
given by 

/ e a/ 3 7 (t • [dpT x <9 7 r]) dS^ = 4irm, (28) 

JdV 

where the integration is performed over the boundary dV of a region V containing 
null-points. 

Conditions (25)- (28) complete the MLR equations in the general case when the 
Jacobian has localized zeroes. 

The above representation involves simultaneous use of Lagrangian variables in 
eqs. (3), (24), (10), (23) and Eulerian ones in (25), making the numerical integra- 
tion of these equations very cumbersome. It is therefore of interest to look for a 
representation formulated in the sole physical space. 

Let us consider the inverse of the mapping a = a(r, t). Using eq. (3), one has 

d t ti + (v ± • V)a = 0. (29) 
From (18), eq. (10) for the magnetic field rewrites 

h = £i jk hQi(a)\Vaj x Va k \. (30) 
Formula (23) for the vorticity in r-variable becomes 

n(r,t) = curl^Vcii) (31) 

where 

V = v (a)+h (a) x W. 
Similarly, the equation of motion (21) for the vector W transforms into 

c> t W + (v ± -V)W = -(j-V)a, (32) 

with initial condition W| t=0 = 0. Here the generalized current j is given by (15). 

These equations are completed by relation (25) and the definition of the normal 
velocity = II v, where the projector II is defined by means of the unit tangent 
vector r = h/h as U a p = 5 a p — r a Tp. They provide a closed system for ideal MHD 
flows, where all the spatial derivatives are taken with respect to r-variables. 
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4 Conservation laws in two dimensions 



The magnetic line representation significantly simplifies in two dimensions where 
the magnetic field lies on the same plane as the flow. It is convenient to introduce, 
instead of the initial position a, the scalar magnetic potential ip defined by 

dip dip 

and a Cartesian coordinate y. 

By fixing ip, we select a magnetic line given by 

dx dy 



dip/dy dip/dx' 

The difference ipi — ip2 is equal to the flux of magnetic field between two lines with 
different values of ip. 

In 2D, ip is a Lagrangian invariant, since it follows from the integration of the 
induction equation (1) that 

^ + (v-V)^ = 0. (33) 

The potential 

iP = iP(x,y,t) (34) 

can then be taken as a Lagrangian marker of the magnetic lines. Solving locally eq. 
(34) in the form y = y(x,ip,t), provides the desired mapping that replaces (2). 

This change of variables, being a mixed Lagrangian-Eulerian one, realizes a 
transformation to a curvilinear system of coordinates movable with magnetic lines. 
In order to implement the transformation from variables (x,y,t) to (x,ip,t) in eqs. 
(33) and (11), we use 

df 1 

-xr = — [ftVi> ~ UVt], (35) 
ot y^ 

df 1 

-7T- = — IfxV^ - UVx], (36) 

df = U_ (37) 
dy 

where derivatives are taken relatively to (x, y, t) in the left hand sides of the above 
equations and to (x,ip,t) in the right hand sides. 

Equation (34) for the magnetic potential then transforms into an equation for 
the magnetic line ifj 

Vt + v x y x = v y . (38) 
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This equation is a kinematic condition. As the equation of motion (3), the dynamics 

of y is prescribed by the velocity component normal to the magnetic field line 

1 



Dt = v±y 1 + Ux where v± = (v • n) and n = (—Vx, !)■ In terms of the new 

V 1 - " 

variables, the magnetic field is given by 



Vx 



X 



h - — h - V- 

which are equivalent to the Cauchy formula (10) for the magnetic field in 2D. The 
derivative in the denominators holds for the Jacobian J. The equation for the 
quantity can be found by differentiating (38) with respect to ip and applying the 
incompressibility condition in the form 

dv x dv x dv y 

This results in a continuity equation for y^, 

d t y^ + d x (v x y^) = 0, (40) 

so that y^ has the meaning of a layer density. 

Another useful relation can be obtained from the equations for the velocity 
components v x and v y that now read 

d t v x + v x d x v x = -d x p + (d f p -j)—, (41) 

y<p 

d t v y + v x d x v y = -(d^p - j)—, (42) 

y<p 

where j = curl h is the current directed along the z direction. It is then convenient 
to introduce 

U = v x + y x v y , 

where y x obeys the equation 

dtyx + v x d x y x + y x d x v x = d x v y 



derived from (38). The function U coincides up to the factor l/yl + y% with the 

1 

velocity component tangent to the magnetic field v T = , U. One easily gets 

dtU + d x (v x U) = -d x (p - v 2 /2), (43) 
that can be viewed as a differential form of the Kelvin theorem. 



9 



Combination of eqs. (40) and (42) gives that w = v y y^ obeys 

d t w + d x (v x w) = -d^p + j. (44) 

To find the analog of (23) in the 2D case, it is convenient to make the change of 
variables y = y(x,ip,t) in the vorticity equation 

d t Q + (v • V)Q = Vj x Vip. (45) 

Substituting relations (35- 37) into (45) and using eq. (38), we get 

Q x j 

d£l + v x d x Q x = — — . 

Vi> 

Equations (40) and (43) provide conservation laws for 2D incompressible MHD. 
They remain valid in the hydrodynamic limit, provided ip is replaced by vorticity 
or by any other Lagrangian invariant. 



5 Possibility of magnetic line breaking 

An important property of the magnetic line representation concerns the compress- 
ibility of the mapping (2), which permits magnetic line breaking. At the breaking 
point, the magnetic field, according to (10), becomes infinite due to the vanishing 
of the Jacobian. As it follows from references [3, 5, 6, 7], the possibility of vortex 
line breaking depends on the space dimension. For two-dimensional flows described 
by the Euler equations, vorticity is perpendicular to the flow plane and therefore 
divv^ = 0. As the consequence, the corresponding mapping is incompressible and 
the Jacobian remains constant. 

For 2D incompressible MHD, the situation is different since the magnetic field 
lies in the flow plane. The velocity can therefore be decomposed into transverse 
and longitudinal components relative to the magnetic field direction. In such a case 
divv^ 7^ and the breaking of magnetic lines is not a priori excluded. 

Let us thus assume that a breaking of magnetic lines occurs. Denote by t = 
t(a) > the positive roots of the equation 

J(a,f) =0, 

and find the minimal value t Q = min a t(a) which defines the first instant of time 
when the Jacobian vanishes. Let a = a be the Lagrangian coordinate of the point 
where this minimum is attained. We first consider that near the singular point, as 
t — > to, the Jacobian behaves as 

J = a(t — t)+ -fijAaiAaj (46) 

where a > 0, 7^ is a positive definite (generically non-degenerated) matrix and 
Aa = a — a . This assumes that the magnetic field does not vanish at the collapse 
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point and in particular that the three vectors dr/dai (i = 1,2,3) lie in the same 
plane, with none of them vanishing. In this case, eq. (10) rewrites 



h = — — — . (47) 

a \to — t) + '-fijAdiAcij 

where b = (h (a) • V a )r| to , ao - This corresponds to a blowup of the magnetic field 

h(a ) like l/(t - t). 

The MHD equations conserve the energy £ given by the sum of the kinetic 

f v 2 r h 2 

Sk = J ~7f~dr and magnetic £h = j energies, that both have to remain finite 

as t — > t • 

Let us the estimate the contribution provided by a possible singularity (47), to 
the magnetic energy 

r b 2 

£ h ~ J j^dr. (48) 

By changing variables from r to a, the contribution to this integral arising from a 
ball of radius R ~ r 1 / 2 where t — to — t and centered in a , rewrites 

gs « b 2 f ^ _ t){ D-2)/2_ (4Q) 

J ar + 'yijdidj 

The size of the retained ball is the largest compatible with the asymptotics. 
The contribution from the other region being most likely finite, we conclude that 
a magnetic field blowup in not excluded in 3D for the assumed expansion of the 
Jacobian. The same conclusion holds if the Jacobian vanishes like (to — t) n at 
the singularity point, with a ball size modified accordingly. At a point where the 
matrix 7 is degenerated with e.g. one eigenvalue Ai being zero, the Jacobian locally 
becomes 

J = a (t -t) + lijafaf + (3a{. (50) 

where a^ holds for the projection of the vector a, transverse to the direction of the 
eigenvector associated with the zero eigenvalue. The contribution of the singularity 
to the magnetic energy then scales like ~ (to — ^) 1 ^ 4 , a behavior which again does 
not contradict the possible existence of a singularity. 

In D — 2, the conclusion can be different. Since the contribution of the selected 
ball to the magnetic energy does not tend to zero as t — > to, a small extension of this 
domain to a ball of size R can lead to a logarithmic divergence ~ B 2 log — — ► 00. 
The divergence becomes more dramatic in the case of a degenerate matrix 7, for 
which £1 ~ (t — t) -1 / 4 . This observation leads us to conjecture that a blowup of the 
magnetic field is probably excluded in two dimensions but not necessary in three 
dimensions. Note that the conservations laws (40) and (43) for the two-dimensional 
problem derived in Section 3, could possibly be useful for a rigorous proof of the 
absence of magnetic blowup. 
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6 Conclusion 



The mechanism for a finite-time singularity addressed in this paper corresponds 
to the breaking of magnetic field lines resulting in a catastrophic amplification of 
the local magnetic field strength. It is worth to notice that this process does not 
contradict the necessary condition for blowup in MHD [11] that represents the 
analog of the Beale-Kato-Majda inequality [12]. According to [11] the velocity and 
magnetic field retain their smoothness on a time interval [0,T] as long as 



Hence a finite-time singularity of any kind must be accompanied by the blow-up 
of fl and Vh. However, this criterion does not exclude a blowup of the magnetic 
field as well. Constraints are nevertheless provided by regularity theorems; a result 
for example states that the solution remains globally smooth if the initial magnetic 
field has a mean component sufficiently large compared to the fluctuations, assumed 
to be localized [13]. This property is a consequence of the fact that only counter- 
propagating Alfven wave packets nonlinearly interact. 

A specific conclusion of this paper is that magnetic field blowup resulting from 
magnetic line breaking is unlikely in two dimensions. Nevertheless, the present for- 
malism cannot capture the behavior near a neutral X-point. Numerical evidence and 
self-similar reductions however indicate that in this case the current amplification 
is exponential in time [14] [15]. 

Furthermore, recent direct numerical simulations of 3D MHD indicate the for- 
mation of quasi-two dimensional current sheets that result in a depletion of the 
nonlinearity strength [16], a mechanism that could prevent singularities. In order 
to validate the blowup scenario discussed in this paper, it is thus of interest to look 
for initial conditions that do not lead to bidimensionalization and has an initial ve- 
locity field whose component transverse to the local magnetic field has a significant 
divergence. 
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